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1.  Introduction 

The  problem  of  determining  the  reliability  of  an  existing  or  proposed 
communication  system  has  received  considerable  attention  in  the  engineering, 
statistical,  and  operations  research  literature  [1 , 5,  9].  For  example,  it  is  important  to 
assess  the  probability  that  a  message  sent  from  a  given  source  arrives  at  its  destination, 
when  the  components  comprising  the  system  are  subject  to  failure.  Unfortunately,  most 
reliability  problems  of  any  substance  are  now  known  to  be  NP-hard  or  #P-complete  [4, 
17, 18,  27].  As  a  result,  researchers  have  focused  on  special  network  structures 
(where  polynomial-time  algorithms  are  possible),  or  have  resorted  to  simulation. 

A  number  of  special  classes  of  undirected  networks  have  recently  been  analyzed 
with  success.  Polynomial-time  algorithms  are  now  available  for  calculating  certain 
reliability  measures  in  series-parallel  [21],  inner-cycle-free  [14],  inner-four-cycle-free  [14], 
and  cube-free  [15]  planar  graphs.  Provan  [16]  has  shown,  however,  that  the  problem  of 
determining  source-to-terminal  reliability  remains  #P-complete  for  the  general  class  of 
planar  graphs.  In  order  to  analyze  more  complex  network  topologies,  the  idea  of  pivotal 
decomposition  [5]  together  with  polygon-to-chain  reductions  [28]  can  be  used  to 
decompose  the  original  problem  into  smaller  subproblems. 

Similar  results  and  tools  are  not  as  available  in  the  case  of  directed  networks.  The 
only  significant  types  of  directed  networks  that  are  known  to  admit  a  polynomial-time 
algorithm  are  "basically-series-parallel"  networks  [2,  3].  Also,  unlike  the  case  for 
undirected  networks,  certain  simplifications  are  not  available  when  carrying  out  pivotal 
decomposition  in  the  directed  case  [1].  Nor  does  there  exist  an  "optimal"  factoring 
algorithm,  such  as  that  demonstrated  for  undirected  networks  [19]. 

This  paper  exploits  the  underlying  algebraic  structure  of  network  reliability 
problems  to  produce  a  general  iterative  algorithm,  applicable  to  arbitrary  directed 
networks.  While  not  polynomially-bounded,  it  is  able  to  generate  reasonable 
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approximations  to  exact  network  reliability  with  a  modest  amount  of  computation. 

2.  Algebraic  Structure 

Suppose  that  G  =  (N,  E)  is  a  directed  network  with  node  set  N  and  edge  set  E,  in 
which  nodes  do  not  fail  but  edges  fail  independently  of  one  another.  The  reliability  of 
edge  e  (the  probability  that  edge  e  functions)  is  denoted  by  pe.  Nodes  s  and  t  designate 
the  specified  source  and  terminal  of  G,  and  we  are  interested  in  calculating  the 
2-terminal  reliability  Rst(G): 

Rst(G)  =  Pr  {there  exists  a  functioning  path  from  s  to  t  in  G}. 

Associate  with  each  edge  k  e  E  a  variable  x^.  Then  the  reliability  polynomial 

Fst(x)  =  Fst(x-| . xr)  associated  with  s  and  t  is  a  polynomial  in  x1 , ....  xr  such  that  if  the 

numerical  values  p-| , ....  pr  are  substituted  for  the  corresponding  variables  x1 , ...,  xr  then 
the  resulting  value  is  the  probability  that  a  functioning  path  exists  from  node  s  to  node  t. 
(If  the  Xj's  were  simply  Boolean  variables,  this  polynomial  would  be  identical  with  the 
structure  function  of  the  system  [5].)  The  reliability  polynomial  can  be  concisely 
expressed  using  two  operations  ©  and  ®  defined  on  polynomials. 

To  begin,  let 

a  ai  ar 

T  -X,  *2  ...xr 

denote  a  monomial  term,  where  each  aj  e  {0,1}.  The  operation  <8>  when  applied  to  terms 
Ta  and  T&  yields  the  term  T°  where  Cj  =  max  {a,,  bj}.  This  operation  is  extended  to 
arbitrary  polynomials  by  distributivity.  The  operation  ©  is  defined  on  polynomials  f(x) 
and  g(x)  using 

f(x)  ©  g(x)  =  f(x)  +  g(x)  -  f(x)  ®  g(x). 

Operations  related  to  ©  and  ®  were  apparently  first  suggested  by  Mine  [13]  and  by 


Kim  et  al.  [10].  More  recently  Gondran  and  Minoux  [8],  and  Shier  [24],  have  formulated 
network  reliability  using  the  operations  0  and  0  defined  above. 

Let  S  denote  the  set  of  all  polynomials  that  can  be  formed  from  monomial  terms  Ta, 
T*3, ...  by  finite  applications  of  the  operations  0  and  0.  Then  it  can  be  demonstrated  [24] 
that  (S,  0,  0)  forms  a  distributive  lattice  with  smallest  element  0  (the  zero  polynomial) 
and  largest  element  1  (the  unit  polynomial).  Suppose  that  Pst  is  the  set  of  simple  paths 
from  node  s  to  node  t  in  G.  Define  the  value  v(P)  of  path  P  to  be  the  product,  with  respect 
to  0,  of  the  edge  variables  along  the  path: 

v(P)  =  ®  II  (xk  :  k  e  p}- 

Then  the  reliability  polynomial  Fst(x)  can  be  expressed  as 

(1)  FstM  =  ®SMP):  Ps  Pst). 

As  an  illustration,  consider  the  standard  bridge  network  in  Fig.  1  having  s  =  1  and  t  =  4. 
Since  there  are  four  simple  paths  extending  from  s  to  t,  equation  (1)  becomes 

F st(x)  =  x^4  0  x1x3x5  0  x2x4x6  0  x2x5  . 

Equation  (1 )  is  just  the  standard  expression  for  the  inclusion-exclusion  formula,  applied 
to  paths  in  the  network  [1].  Expanding  such  an  expression,  using  the  definitions  of  0 
and  0,  and  then  substituting  numerical  values  p^  for  the  corresponding  variables  x^ 
yields  Rst(G). 

A  number  of  techniques  have  been  developed  to  calculate  the  quantity  (1 )  for  a 
general  algebraic  structure  (S,  ©,  0)  satisfying  appropriate  properties  [6,  8].  These 
techniques  can  all  be  viewed  as  different  methods  of  solving  the  system  of  equations 

(2)  z  =  zM0es, 

which  is  linear  in  the  operations  0  and  0.  Here  M  =  (mjj)  is  the  weighted  adjacency 
matrix  for  G,  with  my  =  x^  for  k  =  (i,  j)  e  E  and  my  =  0  otherwise.  Also,  es  denotes  the 
s-th  unit  row  vector.  An  (extremal)  solution  z  to  these  equations  is  known  [6]  to  satisfy 
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Zj  =  FSj(x)  for  all  j  e  N.  Thus  by  solving  such  equations  to  find  Fst(x),  and  hence  Rst(G), 
we  also  obtain  the  (s,  j)-terminal  reliabilities  for  all  j  e  N.  Moreover,  unlike  existing 
methods  for  calculating  2-terminal  network  reliability  based  on  paths  [1 , 1 1],  such 
algebraic  methods  do  not  need  to  first  enumerate  all  simple  paths  joining  the  two 
terminals  [6,  8).  These  paths  are  automatically  generated  in  the  course  of  solving  the  set 
of  equations  (2). 

A  natural  way  of  solving  (2)  is  by  means  of  an  iterative  procedure,  whereby  the 
current  estimate  for  z  is  substituted  into  the  right-hand-side  of  (2),  producing  a  new 
estimate  for  the  solution  vector  z.  In  the  next  section,  we  discuss  a  specific  iterative 
scheme  for  solving  (2)  that  incorporates  special  data  structures  to  streamline  such 
computations. 


3.  An  Iterative  Scheme 

The  basic  idea  of  the  iterative  scheme  presented  here  is  that  of  passing  on,  at  each 
step,  the  information  available  at  node  i  to  each  of  its  neighbors  j,  where  (i,  j)  e  E.  Before 
stating  the  general  iterative  scheme,  the  ideas  will  first  be  illustrated  using  the  network  in 
Fig.  1.  We  will  find  all  reliabilities  Zj  =  FSj(x)  relative  to  the  source  node  s  =  1 . 

In  the  algorithm,  a  polynomial  label  is  associated  with  each  node  j.  At  any  stage, 
l_ABEL(j)  will  be  a  reliability  polynomial  based  on  a  certain  subset  of  paths  from  node  s 
to  node  j.  In  this  sense,  LABEL(j)  corresponds  to  a  current  estimate  of  the  solution  Zj  to 
equation  (2).  Initially,  if  there  is  an  edge  k  =  (s,  j)  e  E  then  LABEL(j)  =  x^.  If  there  is  no 
such  edge  then  LABEL(j)  =  0;  in  the  case  of  the  source  node,  LABEL(s)  =  1 .  Those 
nodes,  apart  from  s,  receiving  a  nonzero  initial  label  are  placed  on  a  list  L.  In  this 
example,  we  have 


LABEL(j): 

L 


1 


0 


xi  x2 

[3,2] 

Now  we  remove  the  "top"  node  i  from  L  and  update  its  neighbors  j  using 
(3)  LABEL(j)  :=  LABEL(j)  ©  [LABEL(i)  ®  xk] . 

where  k  =  (i,  j).  The  above  (i,  j)  update  simply  incorporates  into  LABEL(j)  new  paths  from 
s  to  j  that  use  the  edge  (i,  j).  Any  node  j  whose  label  is  changed  by  (3)  is  placed  on  L  if  it 
does  not  already  appear.  This  step  removes,  in  our  example,  i  =  3  and  updates 
LABEL(2)  =  x-|  +  x2xg  -  x-|X2Xg 
LABEL(4)  =  x2x5 
L  =  [2, 4] 

The  corresponding  network,  with  node  labels  attached,  is  shown  in  Fig.  2. 

At  the  next  step,  node  i  =  2  is  removed  from  the  top  of  L.  Nodes  3  and  4  are  then 
updated,  and  node  3  is  added  to  L: 

LABEL(4)  =  x2x5  +  x-j  x4  -  x^  x2x4x5  +  x2x4x6  -  x2x4x5xg 

-X1X2X4X6  +  X1X2X4X5X6 
LABEL(3)  =  x2  +  X1  x3  *  X1 X2X3  +  X2X3X6  *  X2X3X6 
'  X1X2X3X6  +  X1X2X3X6 
=  x2  +  xix3  -  xix2x3 

L  =  [4,  3) 

This  process  is  continued  until  L  becomes  empty.  At  this  point,  the  polynomial 
label  on  any  node  j  represents  Zj  =  FSj(x).  Table  1  shows  the  final  labels  for  our 
example,  together  with  the  value  obtained  by  substituting  the  common  edge  reliability  p 
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Table  1 


PSj(x) 

Rsj(Q)*  Pk=P  | 

j=1 

1 

1  J 

j  =  2 

X1  +X2X6‘X1X2X6 

p  q  { 

p  +  p^  -  pJ  :• 

j  =  3 

x2  +  x1x3 '  X1X2X3 

p  +  p2  -  p3  l 

j  =  4 

X2X5  *  X1X2X4X6  +  X1 X2X4X5X6 

+  x2x4xs  -  X2X4X5X6  +X1X4 

2p2  +  2p2  -  5p4  +  2p^  1 

: 

-  X1 X2X4X5  +  X1 X3X5  -  X1 X3X4X5 

-  x1x2x3x5  +  xi x2x3x4x5 

V 

1 

-  •: 

The  general  form  of  the  iterative  procedure  is  specified  by  the  following  algorithm, 
where  L  again  represents  the  list  of  nodes  whose  labels  have  been  changed. 


1.  [Initialization] 
for  j*s  do 

If  k  =  (s,  j)  e  E  then  LABELfl)  :=  xk 
else  LABEL(j)  :=  0; 

LABEL(s)  :=  1 ; 

L  :=  [  j:  (s,  j)  e  E]; 
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2.  [Iterative  Step] 

while  L*  [  ]  do 
remove  i  from  L; 
for  k  =  (i,  j)  €  E  do 

T  :=  LABEL(j)  0  [LABEL(i)  <8>  xk]; 
if  T  *  LABEL(j)  then 
LABEL(j)  :=  T; 
if  j  e  L  then  enter  j  into  L. 

Upon  termination  of  the  algorithm,  LABEL(j)  will  be  the  required  reliability 
polynomial  FSj(x).  Notice  that  there  are  several  ways  of  managing  the  list  L.  In  our 
example,  we  treated  L  as  a  queue,  whereby  nodes  are  processed  in  a  FIFO 
(first-in-first-out)  manner.  It  is  also  possible  to  treat  L  as  a  stack,  whereby  nodes  are 
processed  in  a  LIFO  (last-in-first-out)  manner.  The  effect  of  these  two  ways  of  managing 
L  will  be  examined  in  Section  5.  First,  we  discuss  a  number  of  useful  properties  of  this 
iterative  algorithm. 


4.  Properties 

In  this  section  we  make  use  of  the  algebraic  properties  of  (S,  0,  ®)  to  establish 
certain  properties  of  the  iterative  algorithm  presented  in  Section  3.  It  will  be  convenient 
to  denote  the  variable  attached  to  edge  (i,  j)  by  x.  Also,  the  label  on  node  j  at  the  start  of 
step  m  will  be  denoted  by  Lm(j).  Then  the  (i,  j)  update  (3)  of  node  j  after  step  m  is 
expressed  as 

(4)  *-m+i  (j)  =  Lm(j)  0  x  Lm(i). 

Because  the  label  on  node  j  represents  the  sum  with  respect  to  0  of  a  set  of  simple 

i 

I 
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s-j  paths  and  because  this  set  of  paths  can  expand  through  subsequent  updates  (4),  we 
have 


PfQPgrtyJ..  N  k  <  m  then  Lk(j)  ©  Lm(j)  =  Lm(j). 

One  important  simplification  derives  from  the  following  property.  It  states  that  only 
the  "new"  information  N(i)  added  to  the  label  of  i  since  i  was  last  on  L  needs  to  be 
propagated  to  its  neighbors  j. 

Property  2.  Suppose  that  at  step  m  an  (i,  j)  update  is  to  be  performed,  where  the 
labels  on  i  and  j  are  Lm(i)  =  Lk(i)  ©  N{i)  and  Lm(j)  with  k  <  m.  Step  k  represents  the  step 
at  which  an  (i,  j)  update  previously  occurred.  Then  at  step  m+1  the  new  label  assigned 
to  j  will  be  Lm+1 0)  -  1-mG)  ®  x  N(i). 

Proof.  At  step  k,  node  j  receives  the  label  Lk+-|  (j)  =  Lk(j)  ©  x  Lk(i).  Also,  since  k  <  m 
we  have  Lk+1  (j)  ©  Lm(j)  =  Lm(j),  by  Property  1.  Then 

Lm+1  G)  =  ^mO)  ®  x  *-mG) 

=  [LmG)®Lk+1G)]®x[Lk(i)©N(i)] 

=  [LmG)  ©  LkG)  ©  x  Lk(i)]  ©  x  Lk(i)  ©  x  N(i) 

=  [LmG)  ©  LkG)  ©  x  Lk(i)l  ©  x  N(i) 

=  [Lm0)  ®  Lk+i  (j)]  ©  x  N(i) 

=  Lm(j)®xN(i).  ♦ 

Because  the  labels  on  each  node  will  be  maintained  as  fully  expanded 
polynomials  (expressed  using  ordinary  +  and  x),  it  is  desirable  to  know  when  certain  of 
the  terms  in  Lm(i)  do  not  affect  the  label  Lm(j).  The  following  property  provides  one  such 
condition. 

Property  3.  Suppose  that  at  step  m  an  (i,  j)  update  is  to  be  performed  with  Lm(i)  = 

A-|  +  A2  +  ...  +  Av,  Lm(j)  =  B-j  ©  B2  ©  ...  ©  Bw  and  B<|  c  *A<| .  Then  the  updated  label 

Lm+1  G)  =  LmG)  ®  x  \^2  +  +  Ayl- 

Proof.  Let  A  =  A2  +  ...  +  Av  and  B  =  B2  ©  ...  ©  Bw.  Then 


l-m+l(i)  =  f-m(i)  ®  x  l-m(') 

=  [B1  ©  B]  ©  x  [A1  +  A] 

=  [B1  ©  B]  +  [xA-|  +  xA]  -  [xA-j  +  xA]  [B-j  ©  B] 

=  [B"|  ©  B]  +  xA-|  +  xA  -  xAi[Bi  +  B  -  B-jB]  -  xA  [B-j  ©  B] 

=  [Bi  ©  B]  +  xA-|  +  xA  -  xAi  -  xA-j  B  +  xA-j  B  -  xA  [B-|  ©  B] 

=  L m(j)  +  xA  -  xA  Lm(j) 

=  LmG)  ©  x  [A2  +  ...  +  Av] .  ♦ 

Together,  Properties  2  and  3  show  that  certain  "cancellations"  in  the  update  step 
(4)  of  the  iterative  scheme  can  be  predicted  in  advance,  and  thus  unnecessary 
computation  can  be  avoided.  The  next  property  demonstrates  that  the  approximations  to 
RSj(G),  derived  from  successive  labels  at  node  j,  are  monotone  nondecreasing.  The 
notation  RmG)  indicates  the  value  obtained  by  substituting  numerical  values  pr  for  xr  into 
the  polynomial  Lm(j)- 

Property  4.  If  k  <  m  then  R^G)  ^  Rm(j)  • 

Proof.  Since  k  <  m  we  can  express  L^j)  *  Ti  ©  T2  ©  . . .  ©  Tv  and  LmG)  =  ^  © 
T2  ©  . .  .  ©  Tw  ,  where  Tj  is  a  monomial  term  representing  some  path  Pj  from  s  to  j  and 

v  <  w.  Then  RmG)  represents  the  probability  that  at  least  one  path  of  {P-j ,  P2 . Pw}  is 

functioning  and  so  is  at  least  as  large  as  the  probability  R^(j)  that  at  least  one  path  of 
{Pi ,  P2 . Pv)  C  (pi  <  p2-  •  •  -  pw)  is  functioning.  ♦ 

5.  Computational  Results 

Several  examples  will  be  given  in  this  section  to  illustrate  the  efficacy  of  a  version 
of  the  iterative  algorithm  that  makes  use  of  Properties  2  and  3.  The  quality  of  the 
nondecreasing  sequence  of  approximations  to  RSj(G)  will  also  be  examined,  in  particular 
as  this  relates  to  the  discipline  (FIFO,  LIFO)  used  for  managing  the  list  L.  The  iterative 
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algorithm  was  coded  in  FORTRAN  77  and  all  computations  were  performed  using  the 
IBM  3081  computer  at  Clemson  University. 

Example  1 ,  This  network,  having  9  nodes  and  19  edges,  is  taken  from  [20]  and  is 
shown  in  Fig.  3.  There  are  35  s-t  paths  and  5287  noncancelling  terms  in  Fsj(x).  As 
discussed  in  Satyanarayana  and  Prabhakar  [20],  each  noncancelling  term  corresponds 
to  an  "acyclic  subgraph"  of  G.  Despite  its  small  size,  this  example  represents  one  of  the 
most  complex  directed  networks  whose  exact  reliability  has  been  reported  in  the 
literature. 


The  reliabilities  RSj(G)  have  been  calculated  using  our  iterative  procedure  and  the 
FIFO/LIFO  disciplines.  For  ease  of  presentation,  the  reliability  polynomial  Fst(x)  has 
been  evaluated  with  all  p^  =  p  for  the  particular  (s,  t)  pair  indicated  in  Fig.  3;  all  edge 
failures  are  assumed  to  be  independent.  Fig.  4  shows  F(p)  =  Fst(p, ....  p)  plotted  versus  p 
using  the  FIFO  discipline.  As  expected,  the  various  iterations  produce  an  increasing 
sequence  FI , ....  F9  of  reliability  curves  that  converge  to  the  exact  answer  in  9  iterations. 
Each  iteration  produces  a  lower  bound  on  Rst(G)  and  thus  provides  a  conservative 
estimate  for  the  true  network  reliability.  Namely,  the  exact  (s,  t)-reliability  of  the  network 
is  guaranteed  to  be  at  least  as  large  as  the  value  specified  by  the  approximation.  Notice 
that  the  curves  for  the  fifth  through  ninth  iterations  overlap  in  the  figure,  thus  providing 
excellent  approximations  to  Rst(G).  Also  indicated  on  Fig.  4  are  the  cumulative  CPU 
times  (in  seconds)  required  to  complete  the  work  through  the  end  of  the  specified 
iteration.  Thus,  a  total  of  0.636  seconds  were  needed  to  obtain  Rst(G),  whereas  only 
0.061  seconds  were  needed  to  obtain  an  approximation  that  is  virtually  indistinguishable 
over  the  entire  range  0  <  p  <  1 . 

Fig.  5  shows  analogous  information  relative  to  the  LIFO  discipline.  In  this  case, 
twelve  iterations  were  required  before  convergence  was  obtained.  (Several  of  the  curves 
overlap  so  only  10  approximations  are  apparent  in  the  figure.)  Although  the  exact 
answer  was  obtained  in  0.454  seconds  (less  than  the  comparable  time  for  FIFO),  the 


LIFO  discipline  did  not  give  as  useful  a  set  of  approximations  compared  to  the  FIFO 
approach. 

Example  2.  This  network,  with  13  nodes  and  27  edges,  is  derived  from  an  example 
given  by  Martelli  [12];  see  Fig.  6.  It  is  considerably  more  complex  than  Example  1, 
having  70  s-t  paths  and  34,983  noncancelling  terms.  Plots  of  F(p)  versus  p  are  shown  in 
Figs.  7  and  8  for  the  FIFO  and  LIFO  disciplines,  respectively.  Again  it  is  observed  that  the 
LIFO  method  obtains  the  exact  answer  faster  than  the  FIFO  method.  However,  the 
quality  of  approximations  produced  by  FIFO  is  superior  to  those  produced  by  LIFO. 
Indeed,  a  very  close  approximation  to  the  exact  reliability  polynomial  is  obtained  by  FIFO 
in  1 .1 8  seconds,  one-eighth  of  the  time  required  to  find  the  exact  answer  using  FIFO  and 
one-sixth  of  that  required  using  LIFO. 

Finally,  five  random  networks  on  12  nodes  and  30  edges  were  generated  for  test 
purposes.  The  characteristics  of  these  networks,  together  with  the  number  of  iterations 
required  for  convergence,  are  shown  in  Table  2.  In  order  to  compare  the  quality  of  the 
approximations  generated  for  these  examples,  we  have  tabulated  the  CPU  time  (in 
seconds)  required  to  achieve  a  relative  error  of  a%  or  less  (at  p  =  0.5)  in  Table  3.  The 
results  for  Examples  1  and  2  are  also  included. 

In  these  random  examples  the  FIFO  and  LIFO  disciplines  appear  to  be  comparable 
in  terms  of  the  time  required  to  obtain  the  exact  answer.  Again,  however,  the  FIFO 
variant  gives  a  fairly  close  approximation  rather  quickly  and  it  completely  dominates  the 
LIFO  variant  in  this  respect. 
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Table  2 


Network 

#s-t  Paths 

#Noncancelling 

Terms 

iterations 
FIFO  LIFO 

R1 

14 

1,263 

9 

7 

R2 

28 

3,383 

11 

8 

R3 

41 

7,583 

8 

10 

R4 

44 

17,919 

5 

5 

R5 

34 

42,687 

10 

8 

Example 


Discipline 


6.  Conclusions 


This  paper  has  explored  an  algebraic  structure  underlying  certain  network  reliability 
problems.  A  promising  iterative  algorithm  has  been  developed  that  allows  both  exact  and 
approximate  answers  to  be  obtained.  Rather  than  giving  simply  a  single  number,  this 
algorithm  produces  a  reliability  polynomial  that  can  then  be  easily  evaluated  at  any 
particular  input  values  p^ , ....  pr .  Also,  in  the  process  of  determining  Rst(G)  we  also 
generate  RSj(G)  for  all  j  e  N. 

Empirical  results  have  shown  that  the  choice  of  data  structure  (FIFO,  LIFO)  can  have 
a  significant  effect  on  the  relative  efficiency  of  the  procedures  as  well  as  on  the  quality  of 
the  approximations.  Whereas  the  LIFO  approach  frequently  obtains  the  exact  reliability 
polynomial  faster  than  the  FIFO  approach,  the  latter  produces  better  approximations  -- 
ones  that  are  quite  close  to  the  exact  answer  but  are  obtained  in  a  fraction  of  the  time. 

This  desirable  feature  of  the  FIFO  approach  can  be  explained  as  follows,  assuming  that 
the  Pk  are  comparable  in  value.  Under  a  FIFO  discipline,  nodes  are  processed  in  order  of 
increasing  distance  from  s.  Thus,  the  first  time  node  j  is  labelled,  it  is  done  so  relative  to  a 
path  with  the  minimum  number  of  edges.  More  generally,  the  FIFO  approach  ensures  that 
the  "more  probable"  (fewer  edge)  paths  are  incorporated  as  soon  as  possible. 

Subsequent  (longer  and  less  probable)  paths  contribute,  but  not  as  much,  to  the  final 
label  on  node  j.  On  the  other  hand,  a  LIFO  discipline  creates  a  depth-first  rather  than  a 
breadth-first  search  of  the  network,  and  thus  "early"  approximations  can  be  substantially 
improved  by  the  incorporation  of  later  (shorter)  paths. 

The  approximate  solutions  generated  by  the  iterative  algorithm  will  always  produce 
(conservative)  lower  bounds  on  the  exact  solution.  If  greater  accuracy  is  required,  such 
lower  bounds  can  be  used  together  with  a  simulation  approach,  such  as  Fishman’s 
sampling  procedure  [7],  that  makes  explicit  use  of  lower  bounds  to  obtain  improved 
estimates.  Alternatively,  these  lower  bounds  can  be  used  in  conjunction  with  existing 


techniques  that  produce  upper  bounds  on  network  reliability  [22,  25]  to  obtain  an  interval 
that  must  enclose  Rst(G). 

Finally,  it  should  be  emphasized  that  regardless  of  the  list  discipline  used,  some 
relatively  challenging  directed  networks  from  the  literature  can  be  solved  by  our  algorithm 
with  a  modest  amount  of  computation.  In  particular,  one  example  studied  had  70  paths. 
We  are  not  aware  of  any  existing  algorithm  that  has  solved  exactly  a  problem  of  this 
complexity.  While  the  proposed  approach  appears  to  have  potential,  further 
experimentation  will  be  necessary  before  any  firm  conclusions  can  be  drawn  concerning 
its  general  applicability.  In  order  to  solve  larger,  more  realistic  problems  it  may  be 
possible  to  combine  this  approach  with  methods  for  decomposing  the  network  into  more 
manageable  portions  [23,26]. 
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